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Abstract 

A novel model of spontaneous Leptogenesis is investigated, where it takes 
place in the thermal equilibrium due to a background Nambu-Goldstone field 
in motion. In particular, we identify the Nambu-Goldstone field to be the 
Majoron which associates with spontaneous breakdown of (discrete) B — L 
symmetry. In this scenario sufficient lepton number asymmetry is generated 
in primordial thermal bath without having CP-violating out-of-equilibrium 
decay of the heavy right-handed Majorana neutrinos. To obtain the observed 
baryon asymmetry, the neutrino masses are predicted in certain ranges, which 
can be translated into the effective mass of the neutrinoless double beta decay. 



1 Introduction 


Baryon asymmetry of the universe is one of the biggest mysteries in particle physics 
and cosmology. Direct observations show that the universe contains no apprecia¬ 
ble primordial anti-matter, although theories of particle physics treat matter and 
anti-matter in an eqnitable manner [1], To date, the asymmetry between matter 
and anti-matter has been measnred precisely by two independent observations, the 
cosmic microwave background (CMB) measnrement [2] and the measurements of the 
primordial abundances of the light elements [3] predicted by Big Bang Nucleosynthe¬ 
sis (BBN). The baryon asymmetry is often parametrized by the baryon abnndance. 


= 0.0222+°;°“^^ {Planck [2]), 0.0220 ± 0.046 (BBN [3]) , (1) 


or parametrized by the baryon-to-photon nnmber ratio at today’s temperature of 
the universe, 77 ^, 


77b ~ 6 X 10 X 



( 2 ) 


In view of the snccess of the inflationary paradigm, the baryon asymmetry gener¬ 
ated before inflation has been dilnted away, and hence, the baryon asymmetry cannot 
be explained by an initial condition of the universe. As pointed ont by Sakharov [4], 
to generate the asymmetry dynamically, following three conditions are reqnired to 
be fnlhlled in the expanding nniverse; (i) baryon number non-conservation {B vi¬ 
olation), (ii) C and CP violation, (iii) departnre form thermal equilibrium. The 
standard model (SM) could satisfy the all of the conditions, where (i) is achieved by 
quantum effect known as sphaleron [5], (ii) is provided by the CP violating phase in 
the Cabbibo-Kobayashi-Masukawa (CKM) matrix, and the first order electroweak 
phase transition could be realized (iii) if the Higgs boson has an appropriate mass. 
However, the SM fails to generate snfficient baryon asymmetry since the CP violat¬ 
ing phase in the CKM matrix is too small, and the observed Higgs mass 125 GeV 
is too heavy for (iii) [6]. It is therefore of importance to explore the origin of the 
baryon asymmetry in physics beyond the SM. 

Among various possibilities of Baryogengesis associated with physics beyond the 
SM, Leptogenesis is one of the most attractive mechanisms [7]. There, B violation 
is provided by lepton (L) symmetry breaking of the right-handed nentrino mass in 
conjunction with the sphaleron effect, and the condition (ii), the CP violation, is 
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satisfied by the CP-violating phases of the neutrino Yukawa couplings. Then, lepton 
number asymmetry is generated when the right-handed neutrinos in the thermal 
bath decay slowly in an out-of-equilibrium way so that the condition (iii) is satished. 
Leptogenesis is quite attractive since it is naturally achieved along with the seesaw 
mechanism [8, 9, 10, 11, 12] (see also [13]) which explains the tiny neutrino masses 
by the heaviness of the right-handed Majorana neutrinos. 

In this paper, we investigate an alternative model of Leptogenesis which does 
not rely on the out-of-equilibrium decay of the right-handed neutrinos. Instead, we 
consider a model of spontaneous Baryogengesis which has been originally proposed 
by Cohen and Kaplan in Ref. [14, 15], where, in our case, Leptogenesis takes place 
in the thermal equilibrium due to a background Majoron held in motion. In gen¬ 
eral spontaneous Baryogenesis scenario, the background held in motion causes level 
splittings between the matter and the anti-matter via its derivative couplings to 
the current of B symmetry.Then, with efficient baryon violating processes in the 
thermal equilibrium, the non-vanishing level splitting leads to the baryon asymmetry 
of the universe. This mechanism for generating baryon asymmetry works in vari¬ 
ous setups [16, 17, 18], and recently it is discussed that the lepton asymmetry, i.e. 
spontaneous Leptogenesis, can be also achived via anomalous symmetries [19, 20]. 

A notable feature of our model is the use of the Majoron which is the Nambu- 
Goldstone held associated with spontaneous (discrete) B — L symmetry breaking [21, 
22].^^ By remembering that spontaneous B — L breaking is inevitable to obtain 
the mass of the heavy right-handed neutrinos, this choice makes the model more 
economical than the models with additionally introduced held. Moreover, it is also 
advantageous that the derivative couplings of the Majoron to the SM helds are an 
automatic consequence of the spontaneous B — L breaking where they are controlled 
by the B — L charges of the SM helds. 

This paper is organized as follows. In section 2, we explain our setup, and gives 
a formulation to make our notation clear. Our main results are given in section 
3 where numerical analysis of the baryon abundance and some phenomenological 
implications are discussed. Section 4 devotes a summary. 

the literature, the level splitting is often called the effective chemical potential. However, 
this terminology is somewhat misleading, since the chemical potential is defined to characterize the 
thermal bath. In this paper, we use “(dynamical) level splitting” to refer the effective chemical 
potential appearing in spontaneous Baryogenesis. 

^^See, for example, [23, 24, 25, 26, 27, 28, 29, 30, 31, 32] for extensive studies on phenomenological 
and cosmological aspects of light Majoron. 
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2 Spontaneous B — L symmetry breaking and Ma- 
joron oscillation 

In this section we specify the model of our interest, in which right-handed neutrinos 
Nr acquire super heavy Majorana masses via spontaneous (discrete) B — L symmetry 
breaking. For now, let us conhne ourselves to a minimal extension of the SM with 
global U{1)b-l symmetry and one flavor approximation for the lepton sector. To 
achieve the Majorana neutrino mass from spontaneous breaking oi B — L symmetry, 
we introduce a SM singlet complex scalar held, a, with the B — L charge —2. With 
this charge assignment, the Lagrangian is given by, 

C = CsM + \NRi$NR+\d,a\^-[y,NR{L-H) + h.c] 

[oN^Nn + h.c] -V{H,a\ (3) 

V{H,a) = \h{\H\" - vLf + K{\a\^ - 

- vL). (4) 

Here, qn, are dimensionless coupling constants, while Vew and vb-l are 

dimensionful parameters which provide the electroweak and the B — L symmetry 
breaking scales, respectively. The dot products denote the SU{2) invariant products, 
[A ■ B) = A 1 B 2 — A 2 B 1 . Around the B — L breaking vacuum, we parametrize a by 

a = {vb-l + (5) 

where p is a real part of a, and y corresponds to the Nambu-Goldstone boson, i.e. 

the Majoron. After B — L symmetry breaking. Nr acquires the Majorana mass 
Mr = Pn^b-l- Hereafter, we assume pat = 0(1) and omit g^, i.e., vb-l = Mr. 

Since we are mainly interested in the universe where the right-handed neutrinos 
have decoupled from the thermal bath, it is enough to use an effective theory obtained 
by integrating out Nr: 

Ces = (kinetic terms) - Jr - ■ H){L ■ H) + ■ ■ ■ , (6) 

to take the L breaking effects into account. Here, is the fermionic lepton current, 
and rriu = \yy\‘^v‘l.^/M r denotes the mass of the light neutrinos in the SM. To achieve 
this Lagrangian, we have redehned Nr —)■ ]S[^ a^^d L —)■ 

It should be noted that these redehnitions are just a choice of basis in held space. 
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and hence, physical observables do not depend on the choices of the basis.As 
we will see in Section 3.3, the temperatnre Tosc at which the Majoron starts to 
oscillate and spontaneons Leptogenesis takes place snccessfnlly is Tosc ^ 10^^ GeV 
(see Eq. (33)). In the following analysis, we assnme that the right-handed nentrinos 
are heavy enongh, Mr > T^sc ^ 10^^ GeV, so that the effective field theory in Eq. (6) 
is applicable. 

At this stage, Majoron is a massless boson becanse so far it has been treated as an 
exact Nambn-Goldstone boson. It is known, however, that global symmetries might 
not be respected in context of qnantum gravity, and hence violated by gravitational 
effects [33, 34, 35]. Thus, once gravitational effect turns on, Majoron would acquires 
a mass [36, 37] via Planck scale suppressed operators: 

—n 

Om = (n = 5,6,7,...), (7) 

iWpi 

where Mpi ~ 2.4 x 10^® GeV is the reduced Planck mass. We assume that these 
breaking operators appear with 0(1) coefficients which we have omitted here. It 
turns out that the Majoron mass induced by the dimension n operator, \ is 
given by 


m 


(n) 

X 


r\j 


VM^fV 


( 8 ) 


It should be commented that the origin of the explicit breaking terms in Eq. (7) 
may be also attributed to spontaneous breaking of the gauged and hence U{1)b-l 
symmetry broken at around the Planck scale to a discrete B — L {Z 2 n) symmetry. 
In the following section, we discuss how the motion of Majoron leads to spontaneous 
Leptogenesis in this setup. 

For preparation, let us scketch how the Majoron held behaves in the expanding 
universe. We suppose that B — L symmetry is spontaneously broken during/before 
inhation. During inhation, the quantum huctuation of the Majoron is exponentially 
stretched [38, 39, 40], and the Majoron held is settled at some point on its held space. 
After inhation, the Majoron held behaves according to the equation of motion. 


X 3i7x = (9) 

#3lt is also possible to perform computations without these redefinitions. In this case the cur¬ 
rent term does not appear, and instead, the Majoron field appears in front of the dimension five 
operators, ■ H){L ■ H). See the appendix B for details. 
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where H is the Hubble parameter. Here we dehne OqX = X) and take — (m^/2)x^ 
with the Majoron mass m^. In the radiation dominant epoch, the solution of Eq. (9) 
is obtained as 

X{t) = Ji/4(mxt) , (10) 

where, v^Xo is the initial amplitude. This shows that the Majoron starts to oscillate 
coherently when the Hubble parameter decreases down to if ~ m^. During inflation, 
the Majoron is expected to take a field value on its held space, 0 < x/v^ < ttvb-l, 
and hence, the initial value of the motion of the Majoron is xo = 0{Mji)- 

Once the Majoron starts to oscillate, the derivative couplings in Eq. (6) become 
x/{V^Mb)J^, which eventually lead to level splittings between the leptons and the 
anti-leptons as we will discuss shortly. A notable feature of our setup is that the 
Lepton number violating processes automatically present in Eq. (6) as the dimension 
hve neutrino mass operators, and hence, all the necessary ingredients for sponta¬ 
neous Leptogenesis, i.e. derivative couplings to the Majoron and the lepton number 
violating interactions, are equipped with. As we will see in the next section, once 
we specify the neutrino masses, resultant lepton asymmetry depends on only two 
parameters, i.e., the mass and the initial amplitude of the Majoron. In particular, 
certain value of the initial amplitude is automatically chosen due to inflation since 
B — L symmetry breaking occurs before/during inflation in our setup. Although 
the initial value is randomly determined, it is typically to be of the order oi B — L 
violation scale, 0{Mb)- It should be also emphasized that the dynamical level split¬ 
tings violate the CPT-invariance, and hence, baryon asymmetry can be generated 
without satisfying the Sakharov’s conditions exactly. 

3 Spontaneous thermal Leptogenesis 

As we have seen in the previous section, the motion of the background Majoron held 
leads to non-trivial contributions to the kinetic term of the leptons. In this section, 
we hrst discuss the kinematics of leptons in the presence of the background Majoron 
held. Then, we will move on to explore viable parameter regions by solving the 
Boltzmann equation for the lepton asymmetry. 
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3.1 Thermally averaged cross sections in the presence of 
Majoron background in motion 

To see how the level splittings appear, let us look at the kinetic terms of the leptons. 
Here, let us collectively denote the charged leptons and the neutrinos by 'ijj with a 
mass Due to the derivative coupling to the Majoron, the kinetic term of ip is 

deformed by = x/{V2Mji): 

Ckin = ip{i0 - m - n^'y^)ip . (11) 


For Ip having momentum p = {E,'p), the dispersion relation of tp changes to E = 
■E\Jw? + \p\‘^+Px non-zero contribution of /i^, contrary to E^ = 

for = 0. Thus, the term proportional to causes the energy level splittings 
between the leptons and the anti-leptons.^® 

In the presence of the dynamical level splittings, the lepton number asymmetry 
is generated in thermal equilibrium by the L symmetry violating processes via the 
dimension hve operator in Eq. (6). That is, as we will show by solving the Boltzmann 
equations, the L breaking processes smooth out imbalances in the lepton and the 
anti-lepton numbers caused by the level splittings. 

By denoting L = ep)'^ and E[ = {h^, fhe relevant processes for the 

Boltzmann equation of the lepton asymmetry are as follows. 


(a) y-)- e^e^, 

(e) h^hP yy 
{g) h^el yy h~el, 
(k) t-y 


ih) h h yy e^e^, 
(/) yy ejjz/i, 

{h) yy 
(/) h~pL t-y 


(c) yy z/iz/i, 
(i) h^el yy 


(d) yy z/iz/i, 

(j) h~el yy 


We express the scattering amplitudes for the process, for example, h~^h~^ —)■ 
by and —)■ h^h~^ by Ai^- To obtain the Boltzmann equations, it is 

convenient to consider thermally averaged cross sections in which we approximate 
the distribution functions by the Maxwell-Boltzmann distribution, 

(aiv) = f dUidn2dn3dU4e-’^°AT^-E°/TJ2\MY, ( 12 ) 

spins 

#^Here, we discuss the case of the Dirac fermions, although the result can be applied to the 
Majorana fermions. 

^®See the appendix A for detailed discussion. 
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where / = a, a, • • •, and we define dlii = (Ppi/ (27r)^/(2i?°) for the momenta assigned 
by P 1 P 2 —t PsPa in each process.^® Normalization factor is given by = gT^/Ti"^ 
where g is degrees of freedom of corresponding particle. 

The detailed discussion of the cross sections are shown in the appendix A, and 
finally we obtain (a/n) as follows: 

(asv) ~ (don) ^ + 

W-gV) ^ {(Tov) ^ 

where we take massless limit. Other cross sections can be written in terms of them: 

{a^v) = (a^) = = {(^y) = (l/4)(a-n) = (l/4)(a-n), (15) 

Yv) = (a^) = = (a^v) = (l/4)(cT^n) = (l/4)(cT^n), (16) 

(a^) = {cty) = (l/4)(a-n) = (l/4)(a-t;) = (l/4)(a-t;) = (l/4)(a-t;), (17) 
(a^) = (a-n) = (l/4)(a<Tn) = (l/4)(a-t;) = (l/4)(cT-n) = (l/4)(a-n), (18) 

where {aov) ~ ml/{32nvj^). 


3.2 Boltzmann equations 

Before deriving the Boltzmann equations, let us discuss the relations among the 
chemical potentials so that the Boltzmann equations are reduced. First, let us list 
the chemical potentials of the SM particles: 


gauge bosons 

Py, 


Pz, 

Pg: 

matter fermions : 



d'i'Li t 

Pi'Ll 1 



d'Snii 




Pui^i ) 

PuLi ) 

d'd.Li ) 

Pd.Li’’ 


d'UJii I 

f^URi ) 

f^dRij 

Pd-Ri 

Higgs boson : 






where the index i denotes the flavors.In the highly heated thermal bath, the 
gauge bosons have vanishing chemical potentials, and hence, the chemical potentials 

#®Here, we define the chemical potential of the thermal equillibrium so that the Boltzmann factors 
are given by _ ^-^ 12 /'^^ See the appendix A for details. 

#^Here, again we are neglecting electroweak symmetry breaking, and hence, should be 

understood as the ones of the SU{2) x U{1)y gauge bosons, strictly speaking. 
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of the particles and the anit-particles take opposite values with each other. Further, 
we hereafter neglect the flavor mixing for simplicity, so that the chemical potentials 
do not depend on flavors: 


— h'ei) — Ill'Ll 

~ = l^dj^i (dujii = fd'Ujii = fddii- 

The vanishing chemical potential /ivc± = 0 leads to further reductions, = 

hi/i — fisL = f^h+ ~ = 0. For later purpose, we introduce fj,L, hq, and fin which 

denote 


h'fiL — h'U) h'ttL f^d]^ — hQ) h'/l® h'/l+ — jdH- 


( 20 ) 


Neutrality of the universe also puts a constraint on the chemical potentials. In 
general, the total charge of the universe, which is denoted as for a quantum 
number A, is obtained by = Yhi ^diiQAi where Qai is a charge of a particle i, 
and Aui is dehned by An* = — ^.^anti-particie _ 2 ^jT^/ 7 r^(/ij/T) with fXi and 

Qi being the chemical potential of particle i and its degrees of freedom, respectively. 
Here we have again approximated the distributions by the Maxwell-Boltzmann dis¬ 
tribution. The asymmetries between the particle and the anti-particle numbers are 
given by 1/6 and 1/3 instead of 2/7r^ for the Fermi-Dirac and for the Bose-Einstein 
distributions, respectively. Therefore, our numerical analyses are expected to be 
saddled with 0(10)% errors due to the Maxwell-Boltzmann approximation. Thus, 
the total hypercharge, Qy^, is given by 


/^tOt _ 

Vy — 


2T2 


TT^ 




2 + fdd^j 


-\- 2fi 


Ur 




2 + fduL ) 


fdsR 


+ “^(h/i+ + h/i°) 


( 21 ) 


where Ng and Nh are the number of generation and Higgs doublet, respectively. In 
the period when the electroweak interaction is in equilibrium, the Yukawa interac¬ 
tions are also effective in thermal equilibrium^®, which leads to 


fJ'L — — fJ'H — 0, /iQ — + fdH — 0, IdQ — Hdji — fJ-H — 0 , (22) 

where we have used Eq. (20). Altogether, the neutrality condition for the total 
hypercharge gives the relation 

2 T^ 

gy = —[2Ng{fiQ-fiL) + {4Ng + N,,)fiH]=0. (23) 

"^^Strictly speaking, the Yukawa interaction of the electron-type is in thermal equilibrium only 
for T < lO"* GeV. 



Besides, we suppose that the sphaleron interaction is in thermal equilibrium during 
Leptogenesis occurs, which leads to 


-^g(3/iQ + Pl) — 0. 


(24) 


By solving above relations, we can express all the chemical potentials of the SM 
particles in terms of 

_ ^ 4iVg + 3iV/j 

~ 3{4Ng + Nhf^^ ^ 

iNg - Nh l2Ng + Nh 1 

^ 3(4iV3 + iV;,)^^’ ^ ~3iANg + Nhf^^ ^ ~3^^' 

In particular, the baryon number ub = {2T‘^/n‘^)Ng{2^Q + jJ^UR + and the lepton 
number rii = {2T^/n‘^)Ng{2fiB + lien) are given by. 




23Ng + 
3{ANg + Nh) 




(27) 


From this, we also obtain a relation between baryon asymmetry and lepton asym¬ 
metry. 


4(4iV3 + iV0 
23Ng + mh 


(28) 


where ub = —{52/93)nL by substituting Ng = 3 and = 1 in the 

As a result of the above reductions, we are left with only one undetermined 
chemical potential, fii, and hence, we only need to solve the Boltzmann equation 
of fiL- By using the cross sections given in the previous subsection, the Boltzmann 
equation of is given by^^° 


d hl 

dTY 


w 


fJ-L 

Y 


a 


h-x 

T 


(29) 


where w is the wash-out factor dehned by 


Kjapy) Y 
90 sHT ’ 


K ~ 3 X 10^, 


(30) 


ratio of the baryon asymmetry to the lepton asymmetry differs from the ones in the 
literature since the neutrality condition given in Eq. (21) is slightly different due to the Maxwell- 
Boltzmann approximation. 

should be noted that the obtained differential equation is the Boltzmann equation for the 
total lepton charge , and thus, we can safely omit all collision terms other than L violating 
terms: dQ^£^/dt oc {aov). 
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and a is a numerical factor given by a ~ 0.5. Here s and H are the entropy density 
and the Hubble parameter given by 


Jt±T^ i/ = ./hSZl, 

45 ’ V 90 Mpi’ 


27r2 


(31) 


where and g^s are the effective degrees of freedom for the energy and entropy 
densities, respectively. When the temperature is high enough, such as larger than 
0(100) GeV, these two effective degrees of freedom get close to each other around 
g* ~ g*s ~ 100. It should be noted that w is independent of the temperature 
during the radiation dominated period. The Eq. (29) can be further simplihed by 
introducing iii/T = , which can be simply solved by 

(l(T) = - aw I dT (32) 


where denotes the initial temperature to solve the Boltzmann equation. Here, 
Tini takes to be larger than Tosc whose expression is shown in Eq. (33). Since /i^ = 0 
(and thus = 0) at T = our results shown in the next section do not depend on 
the value of After the decoupling of the lepton violating process, i.e. wT 1, 
the lepton asymmetry ends up with /tl/T = ^l(T). 


3.3 Numerical results 

Our goal in this section is to search for viable parameter regions for successful Lep- 
togenesis. After inflation, the Majoron is settled at its initial position, and hence, 
the dynamical level splittings are vanishing in that period, i.e. limt_^ohx(^) 
limt_^o x(t) = 0. Besides, the chemical potentials after inflation are expected to 
be zero, since any asymmetry before inflation has been diluted away by inflation, 
i.e. fiL = fiQ = fiH = 0. With these initial conditions, we solve the Boltzmann 
equation Eq. (32) by the time that the temperature of the universe decreases to the 
sphaleron decoupling temperature Tgph ~ 100 GeV. Around the temperature Tgph, 
the sphaleron rate is sufficiently dumped, and the baryon abundance freezes, which 
is determined by the lepton abundance at that time through ub = —(52/93)ni. 
Notably, free parameters of new physics in Eq. (29) are only and Mb- Thus, by 
remembering that the initial amplitude is of the order of Mb, the baryon abundance 
is given as a function of for given neutrino masses. 

Figure 1 shows the contours of rjB = nB/n,y at today’s temperature as a function 
of and fhu, where ml is dehned by the squared sum of neutrino masses, ml = 
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Figure 1: The baryon-to-photon ratio rjB = nsln^ at today’s temperature is shown as a 
function of and fhi,, where orange lines represent log^g[T b_ 2 ,/Tosc] 

whose temperatures are defined in the text. The left panel of the figure corresponds to 
the case where the initial amplitude of Majoron held is xo = Mr, whereas xo = ttM/j is 
shown in the right panel. 


Here, we are taking the mass diagonal basis of the three neutrinos, so that 
rUyfs appear in the coefficients of the dimension hve operator in Eq. (6). In the hgure, 
we have assumed xo = Mr (left panel) and ttMr (right panel) as typical values. Since 
the wash-out factor is proportional to ml, the net baryon asymmetry is strongly 
washed out when fhj, is large. On the other hand, if is small, sufficient lepton 
asymmetry cannot be achieved since the lepton violating processes are necessary to 
generate the asymmetry from the level splittings. These behaviors of the lepton 
(baryon) asymmetry can be also understood by comparing the temperatures, Tqsc 
and Tr-l, at which the Majoron starts to oscillate, H{Tosc) = lepton 

number violating interactions decouple from thermal bath, H{Tb_r) = ((Jon)T|_^, 
respectively. Suppose the universe is dominated by the radiation, we obtaion 

Tosc ^ 1-7 X , Tr-l ^ 0.33 x gl^'^{(T ov)~^/M r<\, (33) 

where they are typically Tosc ~ 10^^ GeV for ~ 10® GeV and Tr-l ~ 10^^ GeV 
for (cron) 10 GeV The orange lines in Figs. 1 show log^o[7B-L/7ksc]- In the 
region where fhy is large, the lepton number violating interactions strongly couple to 
the thermal bath, and the lepton asymmetry is signihcantly washed out. It turns out 
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that when Tb-l is lower than Tosci h is hard to generate enough lepton asymmetry. 
In the region where fhi, is small, on the other hand, Tb-l is much higher than Tosc-, 
and thus, the lepton number violating interaction is not effective when the Majoron 
held starts to oscillate. As we will see, fhi, is bounded below by experiments. We 
therefore obtain upper limit on fhij for successful Leptogenesis. As a result, we hud 
the allowed neutrino mass ranges for r]B ^ Q x 10“^°, 


< 5.5 X 10“^ eV, 

(XO — Mr), 

(34) 

< 9.1 X 10-2 eV, 

(xo = t^Mr). 

(35) 


It should be noted that the results are saddled with 0(10)% error caused by the 
Maxwell-Boltzmann approximation. 

3.4 Neutrinoless double beta decay 

The allowed ranges of m,y can be translated into the allowed ranges of the neutrino 
spectrum by taking neutrino mass orderings into account, i.e. the normal hierarchy 
(NH) or the inverted hierarchy (IH)’^^^. The NH spectrum corresponds to > 
mu 2 > whereas the IH spectrum to nii.^ > rriy^ > my^. The observed values of 
mixing angles and squared mass differences are given by 

sin2(20i2) = 0.846, sin(20i3) = 9.3 x 10“^ sin(2023) = 1-0, 

rrig = ml^ — ml^ = 7.54 x 10“^ eV^, 

ml = \ml^ -ml\ = 2.47 x 10"^ eV^ (NH), 2.39 x 10“=^ eV^ (IH), (36) 

where we have taken the central values given in Ref. [41]. Then, we have two mass 
spectra as 

NH : my^ = mo, my^ = \Jml + m^, my,^ = \Jrnl + ml, 

IH : my^ = \Jml + ml, m 2 = + m| + ml, my,^ = mo, 

where mo denotes the lightest neutrino mass in each mass ordering. 

Figure 2 shows favored regions of neutrino masses for successful Leptogenesis 
in our scenario, where my-^^, my^, my^^ and ifiy are depicted by three dashed lines 
and black solid line as a function of mo- In the upper two panels of Fig. 2, the 

"^^^Quasi-degenerate spectrum is also possible experimentally. However, such spectrum is not 
favored for successful Leptogenesis in our scenario since wash-out effect becomes too strong. 


(37) 

(38) 
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Figure 2: Possible neutrino mass ranges are shown as a function of mo, the lightest 
neutrino mass. The initial Majoron amplitude is taken to be xo = l^R in the upper two 
panels of the figure, whereas the bottom two panels show the case of xo = The left 

two panels are assumed the NH neutrino spectrum, whereas the IH spectrum is considered 
in the right two panels. The gray regions are disfavored since neutrino masses are too 
heavy/small to produce sufEcient lepton asymmetry. 
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initial Majoron amplitude is taken to be Xo = Mji, whereas the bottom two panels 
represent the cases for xo = The left two panels are assumed the NH neutrino 

spectrum, whereas the IH spectrum is considered in the right two panels. In all the 
panels of Fig. 2, the gray shaded regions depict the disfavored regions for successful 
Leptogenesis, which correspond to the ranges in Eqs. (34) and (35). Therefore, 
the allowed regions exist only in the case that the line of ffij, comes into the white 
region between two gray shed area, which is almost determined by the atmospheric 
neutrino mass scale In the case of Xo = Mr the IH spectrum is rather disfavored 
in all parameter regions, while sufficient lepton asymmetry can be achieved in the 
NH spectrum with the parameter region of mo < 1.3 x 10“^ eV.^^^ On the other 
hand, in the case of xo = both the NH and the IH spectra are consistent 

with successful Leptogenesis where mo < 4.4 x 10“^ eV for the NH spectrum and 
^0 ^ 3-4 X 10“^ eV for the IH spectrum. 

The obtained constraints have implications on the neutrinoless double beta decay 
whose decay width is proportional to a so-called effective mass, |meeP, 

■ (39) 

i 

Here, U denotes the Maki-Nakagawa-Sakata matrix parametrized by 



U 


( I 0 0 \ / Ci3 0 

0 C23 S23 0 1 0 

\ 0 -S23 C23 y y -Sise*^ 0 Ci3 


/ eiai/2 Q Q 

X 0 e*"2/2 0 

^0 0 1 


1 


Cl2 Si2 0 
-Sl2 Ci2 0 
0 0 1 

(40) 


where Sij = sin(6'jj), Cij = cos(0jj), and 5 and ctj are Dirac and Majorana phases, 
respectively. Once we input a neutrino mass spectra, we obtain the effective mass 
as a function of mo as shown in Fig. 3, where the Majorana phases are scanned 
for [0,27r]. It should be noted that though the effective mass also depends on the 
Dirac phase, the dependence is degenerated with those of the Majorana phases. The 
upper limit on mo comes from the CMB observation, ^ 0.9eV [2], which is 

roughly the same for both the NH and the IH cases. The range of effective mass 
|mee| > 0.2 eV has been excluded by the null observations of the neutrinoless double 
beta decay [42, 43]. 

again it should be cautioned that we have 0(10)% errors in our estimations. 
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Figure 3: The effective mass, \mee\, is shown as a function of mg, where the NH and the 
IH cases are depicted by the blue and the green regions, respectively. The upper bound 
for TTT-o comes from the CMB observations, and the upper bound for \mee\ is given by 
Refs. [42, 43]. The red dashed line shows the limit for the NH spectrum in the case of 
Xo = Mr. The red dot-dashed and the solid lines put the upper limits on the NH and the 
IH cases for xo = t^Mr, respectively. 

In Fig. 3, we show the upper limits on mg for the successful Leptogenesis obtained 
in Fig. 2. In the case of Xo = Mr, the IH spectrum is not successful, and hence, only 
the NH case is shown, i.e. mg < 1.3 x 10“^ eV as a red dashed line. In the case of 
Xo = ttMr, the limits, mg > 4.4 x 10“^ eV (NH) and mg > 3.4 x 10“^ eV (IH) are 
depicted by the red dot-dashed and solid lines, respectively. Future measurements 
of neutrinoless double beta decay are expected to reach |mee| = 0(10) meV [44]. 
The hgure shows that if \mee\ > 0(10) meV is conhrmed by the neutrinoless double 
beta decay, almost all the parameter region for the successful Leptogenesis will be 
excluded. 

3.5 Viable models 

Finally let us discuss viable models consistent with cosmological observations. Since 
the Majoron we are interested in is super heavy, m^ > 10^° GeV, Planck suppressed 
operators play important roles not only on their masses as in Eq. (7) but also on its 
decay rate. Here, we consider the following operators as examples: 




( 41 ) 
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to induce a rather large decay rate of the Majoron.^^^ It should be noted that the 
Majoron masses generated by these operators are negligibly small compared to 
in Eq. (7). As discussed below, the power of a in (9^^ is expected to be the same 
as that in Therefore, let us suppose that the Majoron decay is induced by 

(9^^ for example, where Eq. (7) gives to the Majoron mass in this case. On 
the other hand, if dominates the Majoron mass term oX H c:=l Majoron 
would dissipate before the time to oscillate. Let us derive the condition to evade 
such undesired case: 


mix" > x'|h|VM3i~(yVMpi)(TVMpi)~(xVMpi)i7„,,. (42) 


At the time to oscillate, we obtain 


< ~ 2.2 X 10^® 


Mpi \ 

1018 GeV / VlQio GeV 


GeV (43) 


by using Hose = my and x = Xo- We can therefore avoid such situation in most 


cases. 


Since we aim to discuss a connection to a possible ultraviolet completion of our 
model, we constrain ourselves on the cases where both the operators Om and Od are 
induced by a condensate of either a scaler held or a composite held. For example, 
suppose a scalar hied <h with charge q<s> under the gauged U{1)b-l acquires a VEV, 
(<h) = corresponding and are induced from the operators such as 




where we assume > 10, and is a certain high energy scale below Mpp The 
charge is a model dependent parameter, and we demonstrate the cases of = 10 
and 12 in this section, leading to Zio or Z 12 B — L symmetries, respectively. Even¬ 
tually, spontaneous breaking of the resultant discrete symmetry induces Majoron, 
as we have discussed in the previous section. It should be noted that in these cases 
the initial amplitude of Majoron is restricted to smaller range than (0,7r] x vb-l- 
For example, in the Ziq model the remnant discrete symmetry is Z 2 after a ac¬ 
quires a VEV, and hence the Majoron takes an initial held value within the range 
of x/v^ = (0,7r/5] X vb-l- The resultant baryon number asymmetry is, therefore, 
roughly the same as the left panel of the Fig. 1 when xo takes the maximal value. 


#13x^16 decay rate via the derivative couplings is proportional to the neutrino masses, and hence, 
very small. 
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Figure 4: Several cosmological bounds are shown in the plane. The left and 

right panels show the Zio and Z 12 models, respectively. The Majoron never enters the 
thermal bath other than the yellow shaded regions. In the blue regions the energy of the 
Majoron oscillation dominates the energy density of the universe. The deeper blue regions 
represent the range where right-handed neutrino is thermalized after reheating. In the gray 
region Majoron survive until BBN takes place. The black solid lines show the baryonic 
isocurvature fluctuations, and the black dashed lines depict the Majoron masses, , 
given by Eq. (8). 


XO = T^/bVB-L Mr. In this case only NH spectrum is favored as shown in Fig. 2. 
Such restriction is rather relaxed in the Z 12 model. 

Figure 4 shows several constraints on and Mr for two example cases of 
(left panel) and Z 12 (right panel), respectively. In both cases, we require that the 
discrete symmetry is broken before inflation to avoid the domain wall problem. In 
adittion, we also require that the Majoron starts to oscillate well after inflation. 
Thus, the Hubble parameter during inflation, ifinf, is much larger than the Majoron 
mass. In such situation, the quantum fluctuation of Majoron held causes a baryonic 
isocurvature perturbation, which is constrained by CMB observations. In order to be 
consistent with the observations, the Majoron huctuation should satisfy [45, 46, 47]: 

10“® > 5x/X0 ~ i7inf/(27rMij) > Hose/ Mr), wliere Hose = H{Tose) ^ From 

this argument, at least, the condition Hose/{‘^'^M r) ^ 10“^ should be satished. In 
the hgure, we show the value of Hose/ (27rM/j) as the black solid lines. 

In the case of Ziq, the Majoron mass comes from , i.e. in Eq. (8), while 
the x-H-H interaction is induced from (9^^ ~ (M^/Mpj)x|i7p. It is necessary to pay 
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attention whether the Majoron dominates the primordial energy density of radiation 
or not before its decay, since the Majoron decay might dilute the generated lepton 
asymmetry. The temperature at which the Majoron dominates, T^om, depends on 
whether the Majoron is thermalized or not: 

Tdom = — (thermal), w^Tosc (non-thermal). (45) 

TT OiVipj 

The Majoron thermalization occurs via the inverse decay process induced by the y- 
H-H coupling,where the inverse decay rate is given by Tth ~ (^r/^pi)^/ 
which becomes equal to the Majoron decay To aXT ~ The Majoron is therefore 
thermalized when < T^h, where the thermalization temperature Tth is obtained 
by solving the condition H{Tth) = Tth{Tth)- 

The yellow region in the left panel of Fig. 4 satishes < Tth, where the Majoron 
is thermalized by the effect of the decay operator. In that region, however, the 
domination temperature is always smaller than the decay temperature Tq dehned 
by 


Td 


V ^‘g. ) 


(46) 


and thus, the dilution does not occur. In the whole regions other than the yellow 
region, the Majoron does not enter the thermal bath, and the energy density of the 
Majoron oscillation dominates at Tdom as given by Eq. (45). When Tdom becomes 
larger than Td, the Majoron dominates the universe, and causes dilution of baryon 
asymmetry due to entropy production, which is represented by the blue shaded region 
in the hgure. The gray shaded region in the hgure shows that Td is smaller than the 
temperature at which BBN takes place, T ~ 0(1) MeV, where the Majoron decay 
spoils the successful BBN. 

When Tosc > Mr, the right-handed neutrinos are in the thermal bath when 
spontaneous Leptogenesis takes place, which is depicted by the deeper blue shaded 
region in the hgure. In this region, conventional thermal Leptogenesis can take place 
if it is possible, and further lepton asymmetry could be generated even if the dilution 
occurs. Since this possibility is out of our study, however, we do not discuss further. 

Let us hnally mention the Majoron mass rri'x^ represented by the black dashed 
lines in the hgure. In the Ziq model, for example, the Majoron mass is obtained by 

should be noted that the Majoron couplings to the Standard Model particles via the dimen¬ 
sion five operators are highly suppressed by small neutrino masses. 


18 



(5) f5) 

rrix , and the line representing in the left panel of the figure is in the region where 
Majoron dominates the energy density of the universe. This shows that we need 0(1~ 
0.1)% hue-tuning of the coefficient of (9® to obtain viable Majoron parameters for 
a given Mr. The same hgure for Z 12 shows that we need some amount of hne-tuning 
of the coefficient of . 


4 Summary 

We have proposed a new type of spontaneous Leptogenesis caused by the oscillation 
of Majoron held. As exploited by Cohen and Kaplan [14, 15], the mechanism evades 
the Sakharov’s conditions. Therefore, our scenario can work in the case that the 
conventional Leptogenesis does not work, e.g., even when there are not large enough 
CP phases. As a notable feature, all the necessary ingredients for spontaneous Lep¬ 
togenesis are automatically equipped with in our setup, i.e., derivative couplings to 
the Majoron and the lepton number violating interactions. Once the neutrino masses 
are determined, the resultant lepton asymmetry depends on the Majoron mass and 
the initial value of the Majoron amplitude. On the other hand, the dynamical level 
splitting is induced by which is normalized by Mr. Since the initial amplitude of 
the Majoron held is typically to be 0{Mr), is not strongly affected hy B — L sym¬ 
metry breaking scale. Therefore, the sufficient baryon asymmetry can be achieved if 
the Majoron has an appropriate mass. 

To explain the observed baryon number asymmetry, we hnd that the neutrino 
masses are predicted in some ranges since the resultant lepton number asymmetry 
strongly depends on the wash-out effect caused by the dimension hve operators of the 
neutrino masses. As a result, we hnd that spontaneous Leptogenesis rather disfavors 
the degenerate neutrino spectrum. The ehective mass of neutrinoless double beta 
decay is also constrained accordingly. We have also discussed viable models, which 
are Zio and Z 12 models, consistent with cosmological observations. 
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A Cross sections in the presence of the Majoron 
background 

In this appendix we give detailed calculations of scattering cross sections in Majoron 
background. 


A.l Fermions with non-vanishing fi^ 

Let us hrst consider a free Dirac fermion, which corresponds to the limit of the 
vanishing dynamical level splittings, —)■ 0 in Eq. (11). In this case, the Dirac held 
satishes {i$ — m)%lj = 0, and then we expand the solution as 


%Ij{x) = (47) 

^+{x) = j ij-{x) = j , (48) 

where we omit spin indexes. The independent solutions up and vp satisfy the Dirac 
equation for a particle and anti-particle, respectively, and hence, bp and dp are the 
creation and the annihilation operators of the particles and the anti-particles. The 
time dependence of these solutions can be expressed by the energy eigenvalues, 
uJo{p) = a/IpI^ J- as follows: 


where we choose a normalization by recovering spin indexes r and s as 


Of 0 

^P,r^p,s 


= 2Uo{p)Sr 


Of 0 

vJ 

PfT PiS 




(50) 


As a result, we can expand if) in terms of bp and dp^ 


'ip{x) 


d^p 

{27rP2uo{p) 


bpu^pe-^P-^ + div^JP-^ 


p p 


(51) 


where we dehne p = {uo,p). This leads to the momentum conservation of S matrix 
for 2-body process, for example, h'^{pi)h'^{p 2 ) —t e^(p 3 )e^(p 4 ), as 


S = {27rY5{pi + P2 - P3 - Pi) ■ iM, 


(52) 
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where is a scattering amplitude. 

Next let us see the case of non-vanishing where Eq. (52) is slightly deformed. 
In this case, the Lagrangian of a fermion is given by, 

C = ^{i$ — m — (53) 

As discussed in section 2, stems from the Majoron background in motion, and 
thus it depends on time. The situation of our interest is the case at which the 
Majoron starts to oscillate in the expanding universe. The time scale of the Majoron 
oscillation is, therefore, 0(l/m^). The typical time scale of the scattering among 
■^’s is, on the other hand, characterized by 1/T since 'ip is in the thermal bath. 
The temperature at which the Majoron starts to oscillate is roughly obtained by 
H = m^, and hence, 1/T rs-/ I/a/MpittIx is much shorter than the time scale of 
Majoron oscillation. We can therefore treat /i^ as a constant in the cross section 
calculations. 

In the similar way to the previous case, we can take the solution of the Dirac 
equation given by Eq. (53) as 

ijj{x) = ijj+{x) + ip-ix), (54) 

^+{x) = j ^.{x) = j (55) 

It should be noted that these solutions, and vp, have deformed dispersion relations 
given by a;(p) = ±a/|p|^ + nfi -|- /i^ = ±Uo{p} + /i^. The time dependence of the 
solutions is therefore governed by ulp), and hence we obtain 

up= vp= [2uo{p)]-^v^pe-^^^^° . (56) 

It should be noted that and are the one dehned for = 0. This result implies 
that Ai does not depend on /i^, and hence, the S matrix depends on fx^ only through 
the (5-function of the four momentum, i.e. 

S = (27r)^(5'‘(pi +P2-P3-P4 + 2 px) ■ (57) 

where = (/i^, 0) in a certain frame. 

A.2 dependence of phase space integrations 

We now show that the phase space factor does not depend on the effective chemical 
potential at the leading order of First, we again clarify our situation; we are 
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focusing on that the universe is hlled by a background held y which is a function of 
time, and only = {doX, 0) takes place the effective chemical potential. In this case 
the frame taking = 0 is useful for our purpose. However, to calculate the reaction 
rate, i.e., the phase space integration, it is much easier to consider the center of mass 
frame of the initial state. We therefore relate these two frames by Lorentz boosts. 

Suppose the initial particles have momenta pi = {Ei,pi) and p 2 = (i? 2 ,P 2 ), and 
the initial total momentum is dehned by Pjni = Pi + P 2 - In general, when we are 
in the frame where p^ = (p^,0), Pini is not in its center of mass frame. Here, we 
call this frame as the original frame where the momentum set is {Pin^^Px^^}- Each 
momentum is written by 

4?^ = (\/^+|PiniP,Pini), (58) 

pf = (p^,0). (59) 

On the other hand, we can take the frame of = 0, which we call the center of 
mass frame where the momentum set is given by Each momentum 

is written by 


= (Vs.o), (60) 

pP'®’ = ;))|x(\/s+|P,„,| 2 ,-P,„|). (61) 

These two frames, {P/n^^Px^^} and are transformed by the Lorentz 

boost each other: and = A“Vx^\ where 


A = 


/S = 


7 7/5 

7/3^ l + (7-l)S 


A-i = 


7 -7/5 

-7/3^ Il+(7-l)ffi 


P- ■ 

nil 


^+|PiniP 


7 = 


^TlPiniP 


1-|/5|^ 




, (62) 

(63) 


This allows us to obtain the cross section in the original frame after performing 
calculations in the center of mass frame. 

In the center of mass frame we consider the following quantities: 


d^{p2.,Pi) 

I 


d^Ps 


(Pp4 


(27r)32P30(p3) (27r)32P0(p4 

[ d^{p3,P4)'P^\M{pi,P2,P3,pP\‘^, 


{2tt) 6 (Pini - P3 - P4 + 2px), 


spins 


(64) 

(65) 
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where we take massless limit and E^{p) = {i = 3,4). After reparametrizing 

Ps—Px = Ps and PA—Px = Vi) we immediately obtain pa = —ps = —p dne to Pmi = 0. 
The remaining delta fnnction implies i?tot — ~ + Vx) ~ ^ where 

= ^/~S + . By eliminating this delta function, we obtain 


\P\^ = 


EL 


(EL - MP: 






16|^ 


dCMS),2, 


( 66 ) 


n ■ 


where P = \p\n. We finally obtain the expression for 

<^*^(^ 3 ,^ 4 ) as 


d‘^{p3,PA) 


1 _ \p\^dn _ 

4 ( 2 ^)" (£0 + £ 0 ) 1-1 + (£,0 _ E'>)\p^\{h ■ rip 


(67) 


where dhl is the solid angle element of the hnal state momentum. In the original 
frame, there is no specihc three dimensional direction, and hence, the cross section is 
independent of the direction of Pini. We may therefore take Pini = (0,0, p^) without 
loss of generality, which leads \px\ = Px\Pz\/LS and p° = PxVE pD/S. By 
performing the integration over the solid angle, we obtain the cross section. At the 
leading order of /i^, we have 

d^iP3)PA) = d^{p3,pA) + d^{p3,pA) = 7 ^, ( 68 ) 

and hence, the phase space factors do not have Px dependence at the leading order. 
Therefore, we hnd that the Px dependence of the reaction rates, /, should appear 
through the momentum dependence of the squared amplitude after phase space 
integration, 

/~ j d^{p3,PA)^\M{pi,P2)P3)PA)?- (69) 

spins 


It should be noted that in the case where the leptons appear only in the initial 
states, Px dependence would appear through the initial state momenta such as 
■^{.PI)P 2 )P 3 )Pa)- On the other hand, since Px dependence in the delta function 
can be always absorbed by p 3 and pa dependence of the amplitude. Therefore, in 
this case the cross section does not depend on Px since the amplitude is proportional 
to (pi ■ P 2 ), but (p 3 ■ Pa). It is also possible that Px dependence is absorbed by pi 
and p 2 . However, such redehnition of initial state momenta makes the dehnition of 
“actual” lepton chemical potential obscure, and thus, this way is rather complex^^^, 
and we do not employ this manner. 

this way we should introduce the actual lepton chemical potential so that both cases are 
equivalent. 
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B Cross section calculations in other field basis 


We here discuss an alternative and equivalent way to calculate cross sections by 
employing other held basis. When we achieve the effective Lagrangian given by 
Eq. (6), L and is rotated along to y ffat direction, and thus we have the O^x-Jl 
interaction. On the other hand, physical observables should not be responsible to 
this transformation. Therefore, we expect that our result does not change in the 
case that we employ the field basis without the rotation of L. Let us confirm this 
anticipation. 

The effective Lagrangian we are interested in is 




eflf 


(kinetic terms) 




+ ■■• ,(70) 


instead of £efr given by Eq. (6). As mentioned in the appendix A, the time-variation 
of y is much slower than the time scale of the scattering process. This allows us to 
expand y by a reference time as y = yt -|- • • •, and thus we have ~ 

where we can jj,^ as a constant in this calculation.^^® A comparison between £efr, 
leads to a replacement of the couplings from my/{2v‘l^) —)■ without 

the modification of the fermion dispersion relations. This replacement therefore 
modifies the delta function in S matrix, + P 2 — Vz — P 4 ), to, e.g., + P 2 — 

Pz—Pi + 2px) is the case of Eq. (57). Consequently we obtain the same expression 
of the S matrix leading to the equivalent results of the calculation in the appendix 
A. 
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